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Abstract 



We show by means of simple exact manipulations that the thermodynamic per- 
sistent current I(cj), N) in a mesoscopic metal ring threaded by a magnetic flux 
(f> at constant particle number N agrees even beyond linear response with the 
0\ dynamic current Idy(<t>, N) that is denned via the response to a time-dependent 

flux in the limit that the frequency of the flux vanishes. However, it is im- 
d possible to express the disorder average of I,i y ((f>, N) in terms of conventional 

^ Green's functions at flux-independent chemical potential, because the part of 

the dynamic response function that involves two retarded and two advanced 
Green's functions is not negligible. Therefore the dynamics cannot be used to 
map a canonical average onto a more tractable grand canonical one. We also 
calculate the zero frequency limit of the dynamic current at constant chemical 



oo 
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potential beyond linear response and show that it is fundamentally different 
£^ from any thermodynamic derivative. 

*0 PACS numbers: 73.50.Bk, 72.10.Bg, 72.15.Rn 
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I. INTRODUCTION 



The existence of persistent currents in mesoscopic normal metal rings threaded by a mag- 
netic flux (f) has been postulated long time ago [f], although a clear theoretical understanding 
has been achieved only recently [2]. The experimental verification [3]- [5] of this effect has 
given rise to a renaissance of theoretical activity in this field [6]- [11]. Surprisingly, the mag- 
nitude of the measured average current was larger than the available theories could predict. 
The source for this discrepancy between theory and experiment remains controversial. One 
of us has proposed that the long-range nature of the Coulomb-interaction is responsible for 
the large observed current [11], but the issue remains controversial [12] and might eventually 
be settled numerically [13,14]. 

In the present work we shall put aside the fascinating problem of trying to understand 
the Coulomb-interaction, and add some new insights to the dynamic and thermodynamic 
aspects of non-interacting disordered electrons in an Aharonov-Bohm geometry. We consider 
disordered spinless electrons of mass m on a three dimensional thin ring of circumference L 
that is pierced by a flux </>. The system can be described via the stochastic Hamiltonian 



H = H + U(r) , H 
where the momentum operator is 



2m 



P = -V r + — ^e x , <p = j- . (2) 

Here <f> = hc/e is the flux quantum and is a unit vector in azimuthal direction of the 
ring. Throughout this work the charge of the electron will be denoted by — e, and we call the 
azimuthal direction the x-direction. U (r) is a Gaussian random potential with zero average 
and zero-range correlations, 



U(r) = (3) 



*7(r)*7(r') = 7W(r - r') , (4) 

where V is the volume of the ring. We follow here common praxis and assume only short 
range correlations between the random potential. It is possible, however, that an essential 
part of the physics that is responsible for the large observed currents has been lost within 
this approximation, see Ref. [15]. The parameter 7 can be related to the elastic lifetime r by 
calculating the damping of the electrons to lowest order Born approximation. This leads to 
the identification 

7 = 7T~ , 5 
Zttt 

where A is the average level spacing at the Fermi energy. Throughout this work disorder 
averages will be denoted by an overbar. The level spacing A is defined via the relation 

A" 1 = ^ = EKA^ - /^)] , (6) 
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where the particle number N as a function of the dimensionless flux <p and the chemical 
potential ji is for a given realization of the disorder potential U (r) given by 

N(<p,ri = 52f(t a -ri . (7) 

a 

Here f(e a —fi) is the occupation of the exact energy level e„, and /'(e) denotes the derivative of 
the function /(e) with respect to its argument. In a canonical ensemble the function f(e a — ji) 
is not the Fermi function, although for T — > it reduces to a step function irrespective of 
the choice of the ensemble. The energies e a satisfy the single-particle Schrodinger equation 

Hij) a (r) = e a ip a (r) , (8) 

where the ip a (r) are the exact wave-functions. Because H depends parametrically on <p = 
4>/4>o } the energies and wave-functions are functions of <p. Of course, in praxis it is impossible 
to calculate the e a and ip a (r) =< r\a > for fixed and arbitrary potential U(r). Nevertheless, 
the disorder eigenbasis is useful for deriving exact relations between the various physical 
quantities. 



It has been proposed by several authors [7]- [9] that the large measured currents might 
somehow be related to the fact that the rings in the experiment of Levy et al. [3] were not 
attached to external leads. This implies that the number of electrons on each ring is fixed, so 
that for a proper theoretical description one should use a canonical ensemble [7]- [9]. One is 
then faced to the problem of calculating differences between grand canonical and canonical 
thermodynamic averages in mesoscopic systems. Unfortunately, there exist no systematic 
methods to calculate thermodynamic averages at constant particle number, and in praxis 
one has to use some kind of expansion in the fluctuations of the chemical potential [7]- [9]. 
A possible way to circumvent the averaging problem at constant particle number has been 
proposed by Efetov and collaborators [10,16-19]. Let us briefly outline the main ideas of 
this "dynamic approach" to the canonical averaging problem. Suppose that in addition to 
the static flux there is an oscillating component, (f>(t) = (f>o[tp + / du'e~ lw t / n SLp UJ i]. Note 
that in the experiment by Levy et al. [3] the time dependent flux component is given by a 
sine-modulation with amplitude (p ac | and frequency lo ~ 10 — 10 3 i7,2, corresponding to 
Stpu = ^[S(uj + uj ) — S(uj — loo)]. Within linear response theory the time-dependent flux 
gives rise to a time dependent current 8Id y (t,<p, fi) = J dto' e~ lw t ^ n 8I ( iy{to' ^ /i) around the 
ring, with Fourier components given by 

&Idy(v,<P,(J-) = ^-jj- \-^-rj\K{u,(p,n)8(f w + 0(8ipl) , (9) 
L z \mcv J 

where the so-called linear-response function K \u , <p , ji) is given by the Kubo-formula. We 
have rescaled the response function such that the prefactor in Eq.9 is given by 

^(^)=^ = (x)t • A = ^- < u > 

where kp is the Fermi wave- vector and vp = hkp/m is the Fermi velocity. For non-interacting 
electrons K \u , <p , ji) can then be written as 
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K(u,v,ri = K dia (<p, f jL) + K> ara (u>,<p,iJL) , (11) 
where the diamagnetic part is simply given by the negative particle number, 

K dia (i P , l i) = - 1 £f(e a -fi) = -N(ip,fi) , (12) 

a 

and the paramagnetic contribution is 

A^>,^) = E M f{CP ~ ~ /(£ " ~ , (13) 
^ m e a - t fi - u - lO 

with 

P afj = [ dri>* a (r)P x Mr) ■ (14) 



Let us emphasize that Eq.ll is valid in a grand-canonical as well as in a canonical ensemble. 
Of course, in a grand canonical ensemble the electron number in Eq.12 depends on the 
disorder and the flux, while in the canonical ensemble we should solve Eq.7 to obtain the 
chemical potential as a function of <p and N. We shall write //(<£>) for the fluctuating canonical 
chemical potential, and JI for the constant grand canonical one. If an equation is valid for 
both ensembles provided the correct value for the chemical potential is substituted, we shall 
simply use the symbol ji. 

The fundamental assumption of the dynamic approach is [16,17] 



KP ara (0 } (f } JI o ) KP ara (0 } (f 7 ii((f)) . (15) 

This equation tells us that the disorder average of the paramagnetic part of the response 
function is not very sensitive to the choice of the ensemble. Although this replacement is 
claimed to be valid for the paramagnetic part of the linear response function, in a canonical 
ensemble K (0 , <p , /j,(<p)) is the derivative of a flux-periodic function with respect to the flux 
[16,17], so that the average over one period must vanish, 

VW^,/V)) = • (16) 

o 

It follows that at constant N we have exactly 

N= /V^ ara (0,^',/V)) • (17) 
Jo 

Using Eqs.15 and 17 and dividing both sides of Eq.9 by 8<p w , one can write the average 
dynamic current in the limit u — > as follows, 



^-2o 8y w ~ V L J A 



:i8) 



Averaging at constant particle number has now been mapped onto the corresponding aver- 
aging at constant chemical potential, which can be performed by standard Green's function 
techniques. 
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There are some hidden but important subtleties that have been ignored so far. First of all, 
Eq.18 has been derived within linear response theory, which neglects the terms of higher order 
in Stpu in Eq.9. The interpretation of the left-hand side of Eq.18 as an ordinary derivative is 
therefore only justified if the non-linear terms that have been neglected in Eq.9 are consistent 
with the higher order flux-derivatives of the right-hand side of Eq.18. Moreover, even if this 
is the case, it is by no means obvious that the substitution in Eq.15 is indeed correct. At 
least one should be able to estimate the corrections to this replacement. In fact, in Ref. 
[17] Eq.18 was evaluated non-perturbatively in the diffusive regime, with the result that it 
does not precisely agree with the non-perturbative evaluation of an alternative approximate 
expression for I((f } fi((f)) [20]. This discrepancy has never been resolved. In this work we 
shall therefore critically re-examine the validity of Eqs.15 and 18. Our main results are 
(a) that the quadratic response is indeed consistent with the second flux-derivative of the 
equilibrium current at constant particle number, and (b) that the most crucial assumption 
15 is not correct. The result (b) implies that the dynamics cannot be used to map a canonical 
averaging problem onto a grand canonical one, and in praxis one cannot avoid the expansion 
in powers of the fluctuations of the chemical potential [7]- [9], [21]. We derive the correct 
relation between dynamics and thermodynamics at constant particle number and explain 
discrepancies between various approaches that can be found in the literature. Finally, we 
shall discuss the grand canonical linear and quadratic dynamic response functions, which are 
not derivatives of a flux-periodic function. In the case of linear response we re-derive the 
non-perturbative results of Ref. [10] by means of a simple diagrammatic calculation. 

II. THE DYNAMIC CURRENT AND ITS RELATION TO THE CANONICAL 

EQUILIBRIUM CURRENT 

The equilibrium persistent current is given by 



Using the Hellman-Feynman theorem, diagonal matrix elements of the momentum operator 
can be obtained from the flux-derivative of the exact eigen-energies e„, 



Because the three variables /i, N} are related via Eq.7, only two of them are independent, 
and we must specify the variable that is held constant if we differentiate with respect to <p. 
As usual in thermodynamics, the flux derivative of a thermodynamic quantity X at constant 
fi will be denoted by (f^"),u, while the corresponding derivative at constant N is written as 
(fj^)iv- Note that the single-particle levels e a are not thermodynamic quantities, so that the 
flux-derivative in Eq.20 does not depend on whether ji or N is held constant. The following 
manipulations are valid for a given realization of the disorder potential and for a given particle 
number, so that subtleties associated with the averaging procedure [22] are irrelevant. 




a 



(19) 




(20) 
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A. Linear response 



Taking the flux-derivative of Eq.20 at constant /i, we obtain 



dl\ Vfa ( e 2 



L 2 



mc 



V 



Xgc^P,^) 



where the dimensionless grand canonical susceptibility is given by 

L 



2wh 



E 



f(e a - fi) + P aa 



dtp 



dtp 



f22l 



with 



'df(e a -fi)\ , de a 2irh P aa 

dtp ) =f{£a -^ = f{£a -^ — — 



(23) 



Note that we have rescaled the susceptibility such that it can be directly compared with 
the linear response function Kd y (0, tp } /i) defined in Eqs. 11-13. On the other hand, the flux- 
derivative at constant N yields 



'dl\ = V<k 



L 2 \mcV J 

where the dimensionless canonical susceptibility is 



(24) 



T f)P 



(25) 



Note that in a canonical ensemble the derivative N mus t vanish [17]. Comparing 

Eqs. 22 and 25, we conclude that 



Pi 



~ m 



(26) 



Let us now compare Eq.26 with the dynamic response function defined in Eq.ll. In an 
infinite system, the dia- and paramagnetic contributions cancel exactly in the limit u — > 
for any finite disorder. In a mesoscopic system, however, the cancellation is not perfect, so 
that in the static limit the dynamic susceptibility 

Xd y (<P, (*) = lim K(u, tp, fx) (27) 

is finite. In order to see the almost perfect cancellation between the dia- and paramagnetic 
parts of Xdyi^pi l 1 )^ we express Eqs. 11-13 in terms of Green's functions. Defining the advanced 
and retarded Green's functions as usual, 
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we have exactly 



with 



IP „|2 poo Ac r 

A;7>,^,/i) = £^- / —f(e-v){G^e + u)G«(e)-Gi(e)G^e 
H „ m J-oo Ziri 1 ' ' 

afj 



CO 



(29) 



(30) 



and 



IP „ 1 2 poo Ac 

AX ar >,^) = E^^/ T-.Uie-^-fie-u-^G^Gjie-u) . (31) 

a m J-oo ZlTl 

afj 



Using the fact that 



d\a > _ 2ttH . , P (ia 1 
dip ~ L ^ - ' ' V 

dP fyfy 2irh 



p m e a - t fi 



dtp 



L 



1 + 2 J2 



| -Pa/3 | 2 1 



p m c a - cp 



(32) 



(33) 



it is easy to show that in the limit u — > the first and second term in Eq.29 agree exactly 
with the function Xgcif^) given in Eq.22, i.e. 



(34) 



Recall that in a finite disordered metal the energy levels repel each other [23], so that there 
are no degeneracies. As far as A^ ra (u;,<£>, ji) is concerned, we note that the interval of 
integration in Eq.31 vanishes as u — > 0, but the terms with e a = cp in the sum give rise 
to a l/cu-singularity. Cancelling this singularity against the factor of u from the interval of 
integration, we obtain 



afj 



P fj\ 2 

t a ,tp — / ( e a — A*) 



m 



(35) 



where S 6aje/3 is unity if the discrete energy levels e a and cp agree, and vanishes otherwise. We 
conclude that 



ajt/3 



I P I 2 

\ r oifj\ f l ( \ 



afj 



m 



(36) 



Comparing Eqs.26 and 36 and using the fact that S eaj6 = 8 a ^ because there are no degen 



eracies, we conclude that 
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(37) 



In spite of the fact that the above manipulations are very simple and exact, the dynamic 
current calculated in Ref. [17] in the diffusive regime by means of the supersymmetric o- 
model [23] does not agree with the canonical equilibrium current that has been calculated 
in Refs. [7]- [9] via perturbation theory, and by Altland et al. [20] via the non-perturbative 
supersymmetry method. We now explain the origin of this discrepancy. 

It is well-known in the theory of weak localization that disorder averages of products 
of Green's functions of the same type are "harmless" in the sense that they do not involve 
the singular contributions that arise in a perturbative approach due to the famous maxim- 
ally crossed diagrams [24]. With this in mind, the authors of Refs. [10], [16]- [19] (which 
include one of the present authors) have not paid much attention to the disorder average 
Keq ra (0 } <p, fi((f)) } which involves the product of two retarded or two advanced Green's func- 
tions. Note that in the exact disorder-basis 

TSparafn ,„ x f ~ ~ f(tg ~ 



a/3 



\Paf3 12 



-E^^-f'^-^) . (38) 

afj 

The second term in this expression cancels precisely the zero-frequency limit of A^ ra (cu, (f } ji) 
in Eq.35. Of course, if we work at constant chemical potential, then the disorder average of 
Eq.38 can be combined with the diamagnetic contribution K dta (<p, /i) = —N(<p, /i) to yield an 
exponentially small result after averaging, because these terms can be identified with the flux 
derivative of the grand canonical average equilibrium current. In this case the decomposition 
of K para into K'^ ra and K^ 0, is meaningful. However, in a canonical ensemble the first term 
in Eq.38 taken together with the diamagnetic contribution is not exponentially small, so that 
the above decomposition of K' para is not useful. In other words, when evaluating the disorder 
average of Xdy(<P, l^if)), ^ zs n °t allowed to neglect the contribution from — N (<p , /j,(<p)) + 

ATr(o,^M)-^ 

To see how this can explain the discrepancies between the dynamic approach and Refs. 
[7]- [9], [20] consider the Fourier expansion of the equilibrium current I(ip, /i) in Eq.19, which 
is in general of the form 

oo 

I((f,fi) = ^2l n (fi)sm(2irn(f) . (39) 

n=l 

Here the Fourier coefficients I n (fi) are functions of ji and functionals of the disorder potential. 
Similarly, we may expand the relation between N } fi and <p in a Fourier series. Assuming 
that we have solved for ji as a function of <p and N } we have 

oo 

H{V,N) = J2Vn{N)co S {2irn V ) , (40) 

n=0 
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where the ji n depend again on the disorder. Taking the derivative of Eq.39 with respect to 
(f and setting then ji = we obtain 



'dl\ ™ " dIM 

— I = 2^ ln(A i P)) 27TnCOS \ 27Tni P) + _ "a 



s'm(2irn(p) 



(41) 



JV 



Evidently the first term on the right-hand side of this expression corresponds to the term 
Xgc(<P, l^if)) m Eqs.26 and 36, which has been ignored in Refs. [10], [16]- [19]. The crucial 
point is now that both terms in Eq.J f l have the same order of magnitude. Setting //(<£>) = 
~fi + 8fi((f) } where ~fi Q is the disorder average of the zeroth Fourier component in Eq.40, and 
expanding 



dl n (fi) 



d/i 



Sfj, + 0(8^ 



(42) 



fj,=tj, 



we have to leading order in 8fi } 
fdl 



— ] =J2 I n (ii )2irncos(2irn(p) 



N n=l 



+ E 



dl n (fi) 



i dfi 



2irn cos(2irn(p) 



E 



fj,=tj, 

dl n (fi) 



/i - /i + /in' cos(2irn'(p) 



n' = l 



f =1 d/i 



2irn'fj, n i s'm(2irn(p) s'm(2irn'(p) . 



(43) 



Note that the terms with n = n 1 in the double sums contain flux-independent contributions 
that would give rise to an aperiodic current. But these contributions cancel exactly if the 
second and third term are combined, as can be easily seen by writing 



n cos(2irn(p) cos(2irn'(p) — n sin(27rn(,c) sin(27rn'( J c) 

= — [(n — n 1 ) cos(27r(n — n')<p) + (n + n 1 ) cos(27r(n + n')<p)] 



(44) 



Therefore the average over a period of the right-hand side of Eq.44 indeed vanishes, in 
agreement with Eq.16. However, if only the last term in Eq.43 is retained, one obtains an 
aperiodic current that varies linearly with the flux [10]. The second term in 43 explains the 
differences between Ref. [17] and [20]. In the quasi-ballistic regime discussed in Ref. [16] 
the term with ji — /lo in the second line of Eq.43 can be ignored after averaging, because 
disorder averages can be factorized. Moreover, the average current is dominated by the 
diagonal terms n = n' in the double sums. Then it is easy to see that the second term 
is proportional to cos 2 (2irmp) = 7;[cos(4irmp) + 1], while the last term is proportional to 
— s'm 2 (2irmp) = 7;[cos(4irmp) — 1]. Obviously the constant terms cancel and the average 
current due to the last term is exactly half as large as the total current. Unfortunately, the 
factor of two due to the omission of the second term in Eq.43 was compensated in Ref. [16] 
by another mistake, so that the final result was correct and this discrepancy has not been 
noticed. We shall come back to this mistake in Sec. III. In the diffusive regime the averages 
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of products cannot be factorized, so that the term containing fi — fio 1S n °t negligible. Hence, 
even after subtraction of the flux-independent constant [16,17], the current calculated from 
the last term is in the diffusive regime not simply a factor of two smaller than the correct 
result, but its harmonic content is also different. At constant chemical potential the dynamic 
susceptibility is in general not the derivative of a periodic function of the flux. This can be 
seen from Eq.36, since Xgcif^) 1S the derivative of a flux-periodic function, but the term 

|p 12 

J2 a /3 ^t a ,tp l^f f'( e a — fi) is in general not. We shall discuss this term in more detail in 
Sec.IIL 



B. Beyond linear response 



So far, we have shown that in a canonical ensemble the flux-derivative of the thermody- 
namic persistent current agrees with the zero-frequency limit of the linear response kernel, 
i.e. 



lim 



8I dy (u},<p,(j,(ip)) 



'di\ 



(45) 



This suggests that also the higher order flux-derivatives of the equilibrium current at constant 
N agree with the static limit of the corresponding higher order dynamic response functions. 
We now show explicitly this agreement for the quadratic response. Writing 



hy = SI dy + S 2 hy + 0(8(f 3 ) 

the second order contribution to the dynamic current is 

8 2 I dy (t^,fi) = l - {^pj ^ j du l du 2 K( 2 \u u u 2 ^^)e-^ + ^ t IH^ 1 8^ 
with the quadratic response kernel given by 



(46) 



(47) 



K {2) {tO U L0 2 ,ip,fl) 

~ [IT 




Pqp8p a f(ep - fi) - f(e a - jl) 
m C{3 - e a + ui + lO 

Pqp8p a f(ep - fi) - f(e a - fi) 
m - e a + ui + u 2 + iO 



L- 7 



e a + ui + u 2 + iO 



f(tp - ft)- f(e a - ft) /(e 7 - ft)- f(tp - fi) 



ep - e a + cut + iO 



e 1 - ep + uj 2 + iO 



(48) 



This expression can be derived for example by non-equilibrium Green's function methods. 
It gives the stationary state after adiabatically switching on the periodic time-dependent flux 
components. The first and second term contain one current-vertex and one charge- vertex 
and vanish exactly. The first term accounts for the first order variation of the total number 
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of electrons which, of course, is zero. The second term gives the variation of the electron 
distribution functions due to the term in the Hamiltonian which is quadratic in the time- 
dependent vector potential. As this term is proportional to the operator for the total number 
of electrons, it commutes with the Hamiltonian and gives no contribution in this dynamic 
problem. The quadratic response kernel is therefore given by the third term containing three 
current-vertices, which can be written as 

tA^)(, , , , ,„ ,,\ o f^h\ x . Pgf j Pf jl P 1 a \ f{e a - fi) 

\ P J Tpq m {{e a - ep - cut - iO)(e a - e 1 - cut - uj 2 - i0) 

+ f(tp-v) 

(e/3 - e a +u 1 + iO)(tfs - e 7 - u 2 - iO) 

+ fih^fl ) . (49 ) 

(e 7 - e a + ui + u 2 + «0)(e 7 - tp + u 2 + iO) J 

Defining the dynamic second-order susceptibility as the static limit of the second-order re- 
sponse kernel, i. e. 

XdJ&'f 1 ) = % A " (2) (^i^2,</^) , (50) 

w 2 -S-0 

we obtain from Eq.49 

\ L J ap 1 771 



e l3 



^ e ~* ^ i ol = (3 7^ 7 and cyclic permutations 

, a = (3 = 7 

Rearranging by cyclic permutations we finally obtain 

(2)/ x r- ^UX v^/P«flPfl 7 P 7 a J \ ^a. ) 

Xlly'iV^) = "6 T L 



\ "* -P«a|-P«7| /( e 7 Ap fi^a Ap I 

J ' ( J 

The prime at the first sum indicates that a ^ (3 ^ 7 7^ a. 

Now this result is compared with the canonical equilibrium second-order susceptibility 
defined by 

/ d 2 l\ evp 2ir c 9W . . 

M w = -rT^") • (5:i) 



In a canonical ensemble 



'd 2 l\ -e d 2 P aa 

-7T-;) = — rE a r /( e » ~ AO > ( 54 ) 
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so that 



(55) 



2nh^ dtp* 

From second order time-independent perturbation theory for the non-degenerate case we have 

d 2 



IP I 2 

l i a l 



L ) m 2 y (e a — e 7 ) 2 



\a > 



(56) 



+2 £ 



E 



P(3a. Paa 



1/5 > 



<9 P«a g / 271"/?. \ 1 J v Pa[3 P[3~/ P'ya Paa\Pa-y\ 

d( f 2 \ L ) m 2 \ p(¥a) (e a - ep)(e a - e 7 ) (e a - e 7 ) 2 

We conclude that the second order canonical equilibrium susceptibility is given by 

(2)/ \ r- / 27t/?.\ 1 ^— v I ^— v P<y.B Pp-y P-ya ^"aa: | -fary | I 

V L J m a p&a) l e « - e /?)l e « - e 7) l{ -L a ) l e « ~ e 7) 

(58) 

This expression is easily seen to agree with the right-hand side of Eq. 52. Thus 

x£W) = x£ 2 W) > (59) 

i. e. the equality of the canonical equilibrium susceptibility and the dynamic susceptibility 
holds in linear and quadratic order. We suspect that this agreement holds for all higher 
derivatives. We would like to emphasize, however, that we have not proven this agreement 
to all orders in perturbation theory. 




III. DYNAMICS AT CONSTANT CHEMICAL POTENTIAL 

In this section we shall study the average dynamics of mesoscopic rings at constant 
chemical potential. Experimentally this corresponds to the dynamic response of an ensemble 
of rings that are somehow coupled to an external reservoir, so that the number of electrons 
can change as the flux is varied. Although it is hard to imagine how such a situation can be 
realized experimentally, the following calculation is instructive from the theoretical point of 
view, because it shows that the differences between canonical and grand canonical ensembles 
in mesoscopic rings manifest themselves in the dynamics perhaps even more drastically than 
in equilibrium properties. 
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A. Linear response 



We first discuss the linear response within a simple perturbative approach. Combining 
Eqs.34 and 36, setting ji = Ep = const, and using f'(e a — Ep) — > —S(e a — Ep) as T — > 0, we 
obtain in the zero-temperature limit for the average grand canonical dynamic susceptibility 
defined in Eq.27 



Xdy (<P, E F ) = -N(<p, Ep) + ATr (0, E F ) - £ S ea ^^^S(e a - Ep) . (60) 

afj 

Let us now examine this expression in the diffusive regime, where the elastic mean free path 
I is small compared with the circumference L of the ring, but the localization length £ Ml 
is still larger than L. Here M = ^ kp ^^ is the number of transverse channels, where L±_ is 
the transverse thickness of the ring. Due to an almost perfect cancellation between the dia- 
and paramagnetic contributions, the sum of the first two terms in Eq.60 is of order e _i/ ^, i.e. 
exponentially small in the diffusive regime, so that only the last term survives. In Ref. [10] 
this term has been studied by means of the non-perturbative supersymmetry method, with 
the result that its average over the flux does not vanish. In subsequent work [16,17] this flux 
average was subtracted again, and not much attention was paid to this term. Because the 
machinery of the supersymmetry- method is physically not very transparent, let us give here 
a simple diagrammatic derivation of this contribution. 

The last term on the right-hand side of Eq.60 involves a product of two matrix elements, 
but only a single energy denominator. Because the usual Green's functions in the momentum 
or real space basis combine always one matrix element with one energy denominator, such an 
expression cannot be directly written in terms of Green's functions. (Note that 8(e a — Ep) = 
~^i[Ga(Ep) — G^(Ep)], so that a Dirac-$ should be counted as an energy denominator.) We 
therefore follow Ref. [16] and smooth out the Kronecker-$ by replacing 

J2Se a ,e,\P a p\ 2 S(e a - E F ) -> A ^2 P aa Ppp8(e a - E F )S(e p - Ep) , (61) 

where the Jacobian for replacing the Kronecker-$ by the Dirac-$ is simply given by the average 
level spacing, see Eq.6. Note that for spinless electrons in three dimensions A = ir Ep / '(AM), 
where A = kpL. The Jacobian used in Eq.3.11 of Ref. [16] was too large by a factor of two, 
which led to an exact compensation of the mistake due to the omission of the second term 
in Eq.43, and the final result was correct. 

Given Eq.61, the number of matrix elements matches again the number of the energy 
denominators, so that we may use standard Green's function techniques to average this 
expression. Going to the momentum basis |k >, we have to calculate 

Xdy (^E F ) = -A (^)\ h -^[Gi^Ep) - G^(Ep)] {G£, y (Ep) - G^ y (E F )} , 

(62) 

where k x = k x + and 
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Gi£(t) = EMVrMGy R (e) (63) 

a 

is the Green's function for a given realization of the disorder in the momentum basis. The 
wave-functions are V'a(k) =< k|a >. Because we are now working at fixed chemical potential, 
only the combination involving the product of an advanced and a retarded Green's function 
has to be retained in the diffusive regime. The dominant diagrams that determine the disorder 
average in Eq.62 are shown in Fig.l. These diagrams correspond to the following expression, 

X*y(<P,E F ) = -J4iE [Gi(E F )} 2 {G*(E F )} 2 [C(k + k', 0) + D(k - k', 0)] , 



kk' 



(64) 



where the averaged Green's functions are 



G i(c)= f _nA_-n_ > G ^(e)= ^ , (65) 

2m 2r 2m 2r 

and the Cooperon- and diffuson propagators are 

C(w) = i X . ■ (66) 

hT>(q + -j^e x y - iuj 

D(q,u,) = i \ . . (67) 

Here T> = v F £/3 is the diffusion coefficient, and the parameter 7 is given in Eq.5. The 
diffuson contribution in Eq.64 has already been written down in Ref. [16], although it was 
not evaluated there. We would like to emphasize that the diffuson appears in Eq.64 on equal 
footing with the Cooperon. In contrast, the leading weak-localization correction to the Drude- 
formula involves only the Cooperon [24]. The reason for this difference is that the diagrams 
in Fig.l are not conventional conductivity loops, but resemble more Hartree diagrams, which 
have the unusual feature that there is no energy integration associated with the Hartree loops. 
This is due to the fact that Eq.64 involves only propagators at the Fermi energy. From Eqs.66 
and 67 it is clear that D(0,0) and for <p = also C(0,0) are formally infinite. Within the 
framework of perturbation theory this singularity can be cured if contributions from higher 
order diagrams involving more diffusons and Cooperons are resummed [16]. Alternatively, 
we may simply use the results of the non-perturbative supersymmetric cr-model calculation 
[23], which imply that the correct way to regularize the singularity is to shift the frequency 
according to [25] 

uj^uj + iT , T = — . (68) 

7T 

It is important to stress that terms with more than one Cooperon and diffuson are only im- 
portant as far as the zero-mode is concerned. For all other modes the perturbative expansion 
in powers of Cooperons and diffusons is controlled by the small parameter where E c = jj 
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oc so that this expansion is good as long as 



is the Thouless energy. Note that jjh = ^-jfj ^ ^ 
the size of the system is small compared with the localization length £. In this regime it is 
sufficient to retain only the leading term involving one Cooperon and diffuson. 

Using the regularization 68 and assuming that Lj < 
forward. The final result is 



j_l ro £, the evaluation of Eq.64 is straight - 



Xdy(<P,E F ) 

where we have defined 



A 



(2ir) 3 M 



E 



i 



i 



n 2 + 7 (n + 2(f) 2 + 7 



A 



7 



(69) 



(70) 



The series in Eq.69 can be summed exactly. Using Eqs.9 and 10, we finally obtain for the 
dynamic current in the static limit 



lim 



8I dy (u;,cp,E F ) 



( evp 



1 



7T 



'1 + e" 2 ^ N 



V L J (2tt) 2 M v^y V 1 - e~ 2v ^ / 



'I 



-27rv^\2 



1 — 2 cos(47r(,c)e' 



+ e 



(71) 



Expanding this expression for small 7 and using 7 1 = |(27r) 2 My-, we have to leading order 



lim 



8I dy (u;,cp,E F ) 



8tpu 



ev F 



2£_ 
~3L 



1 



A 



A + 2ttE c [1 - cos(4tt(^)] 



(72) 



Note that the first term in the square brace is due to the diffuson pole, while the second 
term is due to the Cooperon. At zero flux both terms cancel, so that at <p = the linear 
response function vanishes. For |27r(,cmod7r| ^> \J^~ the Cooperon contribution is completely 
negligible. Writing Eq.72 in terms of the diffusion coefficient T> = v F £/3, we obtain in this 
regime 



lim 



8I dy (u } (f } E I 



8<p u 



2eV 



|27r(,cmod7r| ^> 



A. 



(73) 



This expression agrees exactly with Eq.18 of Ref. [10] (taking into account that we are 
working here with spinless electrons). Note that Efetov [10] has obtained this result by means 
of the non-perturbative supersymmetry- method, which is in principle exact. Therefore the 
approximations in Eqs.61 and 68 are justified a posteriori. In Ref. [17] the right hand side of 
Eq.72 is integrated over <p after the average over the flux has been subtracted, and the result 
is interpreted as the canonical equilibrium current. From the analysis of Sec. II it is clear that 
such an interpretation is not correct. From our derivation it is also obvious that the value 
of the response function in the regime |27r(,cmod7r| \J~^- is completely determined by the 
pole of the diffuson. Because the existence of this pole is a consequence of particle number 
conservation, Eq.72 should be not very sensitive to inelastic processes [10]. Note, however, 
that the zero-mode which is responsible for the large value of the linear term exists only in 
an isolated ring. If the current is measured via leads that are attached to the ring such that 
the azimuthal symmetry is broken, the zero mode should be omitted. 
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B. Beyond linear response 



We now calculate the leading non-linear correction to Eq.71. Writing Eq.48 in terms of 
Green's functions, we obtain 



K^\u u u 2 ^E F ) = -A f^)^ W79 r A 

TJl y Li J — oo ZTTl 



E F ) - f(e + u 2 - E F )]G R (e + u x + u 2 )G R (e + u 2 )G A (e) 
+ [f(e + u 2 - E F ) - f(e + u 1 +u 2 - E F )]G R (e + + u 2 )G A (e + u 2 )G A (e) 
+f(e - E F ) [G A (e)G A (e - ^)G A {e - - u 2 ) 

-G R (e + u l+ u 2 )G R (e + u 2 )G R (e)\} . (74) 

In a grand canonical ensemble we know that disorder averages involving the combinations 
G A G A G A and G R G R G R are exponentially small in the diffusive regime, so that only the first 
two terms in Eq.74 survive after averaging. We emphasize again that in the canonical en- 
semble this is not the case. The response function K^ 2 \— cu, cu, ip, E F ) in the frequency regime 
A ^ io ^ E c has been studied in Ref. [26]. Here we are interested in the limit Ll>i } Ll> 2 — >■ 0, 
where we encounter the same problem as in the case of linear response: the number of matrix 
elements exceeds the number of energy denominators by one. In order to write this expres- 
sion in terms of Green's functions, we introduce again an additional energy denominator 
via a smoothing procedure similar to the one used in Eq.61. In this way we obtain in the 
momentum basis 

2tt 1 



kk'k" 



X \pY^F) ~ G k,k( EF )} l G k',k"( EF ) G V',k'( EF ) ~ G k' ,k"{ E F) G k" ,k'{ E F)\ ■ ( 75 ) 

The prefactor of 3 is due to the three cyclic permutations in Eq.51. On the other hand, if 
we differentiate the average linear response function Xdyi^pi E f) in Eq.62 with respect to (f } 
using the exact identity 

^~ G ky( e ) = ~~T X! G ky'( e )K G k"y( e ) > ( 76 ) 
dtp m Li k „ 

we obtain an expression almost identical to Eq.75, the difference being that the prefactor of 
3 is replaced by a prefactor of 2. We conclude that in the diffusive regime 

3 d 



X^,E F ) = ~Xd y (v>,E F ) . (77) 

Note that the factor of | does not appear in the corresponding equation for the canonical 
ensemble, see Eq.59. According to Eq.77 the large coefficient of the linear response term in 
Eq.73 does not appear in the quadratic response, so that outside the narrow regions where 
2(f is close to an integer the leading correction to the linear response is small, even if Stp^ is 
of the order of unity. 
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IV. CONCLUSIONS 



In this work we have clarified the relation between the dynamic current and the thermody- 
namic current at constant particle number, and have settled a controversy that has remained 
unresolved for the past three years. We have shown that the fundamental assumption of 
the dynamic approach, Eq.15, is not justified, because the flux dependence of the chemical 
potential in the disorder average K£q ra (0, <p, /^(f)) cannot be neglected. Although this part 
of the response function involves only disorder averages of the form G A G A and G R G R } its 
contribution to the average canonical persistent current is of the same order of magnitude as 
the contribution from the G^G^-term that is usually retained. Thus, one of the most basic 
properties of disorder averages of products of Green's functions in the diffusive regime does 
not apply to mesoscopic systems at constant particle number. Hence, also in the calculation 
of dynamic properties one cannot avoid the problem of expanding in powers of the fluctu- 
ations of the chemical potential, assuming that such an expansion is allowed. In the light of 
this result previous calculations should be critically re-examined. 

We have also calculated the leading corrections to the linear response functions and have 
shown that in a canonical ensemble the zero-frequency limit of the quadratic response can 
be obtained from the flux-derivative of the corresponding linear response function. In a 
grand-canonical ensemble, however, this is not the case. Finally, we have shown that non- 
perturbative results obtained via the supersymmetric cr-model can be exactly reproduced by 
combining diagrammatic perturbation theory with a simple regularization prescription of the 
Cooperon and diffuson propagators. 
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FIG. 1. Dominant diagrams that determine the average dynamic susceptibility Xdy{<-Pi Ef), see 
Eq.64. A solid arrow with label R or A denotes an averaged retarded or advanced Green's function 
with energy fixed at Ep. The vertex in (a) represents the Cooperon, and the vertex in (b) represents 
the diffuson. The small black circles denote current vertices. It is understood that there is no energy 
integration associated with the loops. 
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